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Elastic turbulence in shear banding wormlike micelles 
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We study the dynamics of the Taylor-Couette flow of shear banding wormlike micelles. We focus 
on the high shear rate branch of the flow curve and show that for sufficiently high Weissenberg 
numbers, this branch becomes unstable. This instability is strongly sub-critical and is associated 
with a shear stress jump. We find that this increase of the flow resistance is related to the nucleation 
of turbulence. The flow pattern shows similarities with the elastic turbulence, so far only observed 
for polymer solutions. The unstable character of this branch led us to propose a scenario that 
could account for the recent observations of Taylor-like vortices during the shear banding flow of 
wormlike micelles. 
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Wormlike micelles are elongated, polymer-like struc- 
tures resulting from the self-assembly of amphiphilic 
molecules in aqueous solution P, 0j- In contrast to 
regular polymers, they continuously break and fuse, 
providing an additional relaxation mechanism. In the 
fast breaking regime, Cates' reptation-reaction model Q 
predicts that wormlike micelles solutions relax mono- 
cxponentially with a single time tr ~ ^/TbT r , where r& 
and r r are the breaking and rcptation times. 
In addition to their structural analogy, polymers and 
wormlike miccllar solutions can exhibit flow instabilities 
when submitted to even moderate shear rates. In par- 
ticular, many wormlike micelles solutions have been ob- 
served to undergo a shear banding transition. Under 
simple shear, the base scenario is the following [l], 0] : 
below a critical shear rate 7; ~ 1/tr, the flow is ho- 
mogeneous. Above 7j, the system becomes mechanically 
unstable. A phase of lower viscosity nucleates, inducing 
a banded state in which the initial viscous phase and the 
fluid phase coexist at constant stress. Increasing the im- 
posed shear rate only affects the relative proportions of 
each bands, up to a second critical value jh, where the 
high shear rate phase entirely fills the flow geometry. Be- 
yond 7/j, the homogeneity of the flow is recovered. Albeit 
extremely well documented, the shear banding instabil- 
ity is still not fully understood 0-0] • In particular, it 
has been shown recently in cylindrical Couette geometry 
that contrary to the usual view, the shear banding flow 
may not be purely one-dimensional, but instead is orga- 
nized into Taylor-like vortices stacked along the vorticity 
direction. These cellular structures are mainly localized 
in the high shear rate band and exhibit a complex dy- 
namics depending on the applied shear rate . 
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Recent theoretical developments have tried to rationalize 
the 3D flow, by involving an interfacial mechanism driven 
by a jump in normal stresses across the interface between 
bands [f| 0. Besides, complex 3D flow structures are 
also well-known to develop at low Reynolds number Re 
in polymer solutions [10h12| ]. The underlying instabili- 
ties are a consequence of the non-Newtonian stress field 
that is created in the flow due to elasticity. The im- 
portance of the clastic nonlinearity is expressed by the 
Weissenberg number, defined as Wi = 775? = N\/a, 
where N% is the first normal stress difference and a is 
the shear stress. To a good extent, those two definitions 
are equivalent (see Supplementary Figure). Thereafter 
we only use Wi = jtr. Wi acts as a control parameter 
analogous to Re in Newtonian fluids. When increasing 
Wi, a viscoelastic fluid is likely to undergo a transition 
from ID flow to various coherent 3D flows, which would 
eventually lead to turbulence [Rj - floT ]. For polymer so- 
lutions flowing in curved geometries, elastic instabilities 
are triggered above a threshold that follows a general 
criterion established by Pakdel and McKinley and corre- 
sponding to values of Wi ranging typically between 1 and 
10 [l(| Ha]. In wormlike micelles, the onset of the shear 
banding regime is characterized by Wi ~ 1 [l[. Hence, 
on the stress plateau and above, the Pakdel-McKinlcy 
condition is satisfied, suggesting that the elastic instabil- 
ity could be responsible for the 3D shear banding flow [1] . 
In this letter, we investigate the flow behavior of the high 
shear rate branch (Wi > Wih = 7ft Tr) of a shear band- 
ing wormlike micelles system in Taylor-Couette geome- 
try. We show that once the induced phase fills the whole 
sample (Wi > Wih), the vortex structure is destroyed 
and the ID structure of the flow is recovered. Never- 
theless, above a critical threshold Wi c , we observe the 
growth of another instability leading to a random flow 
state, that presents all the features of the elastic turbu- 
lence. This indicates that the induced phase in shear 
banding wormlike systems can indeed be subject to elas- 
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FIG. 1: Normalized shear stress a/Go as a function of the 
Weissenberg number Wi = ^tr. The data are gathered for 
increasing Weissenberg numbers from to 34. Inset : Views 
of the velocity gradient-vorticity plane for different Wi. The 
inner and outer cylinders are respectively on the left and right 
sides of the pictures, (a) Wi < Wit : homogeneous flow of the 
low shear rate phase, (b) Wii < Wi < Win ■ banded state. 
The interface between bands presents oscillations that scale 
with the size of the high shear rate band, in which vortices 
stacked in the vorticity direction develop, (c) Wih < Wi < 
Wi c : homogeneous flow of the high shear rate phase, (d) 
Wi > Wi c : turbulent flow of the high shear rate phase. 



tic instability. From this observation, we propose a sce- 
nario that points to a bulk elastic instability as the un- 
derlying mechanism for the vortex flow in the banding 
regime. 

The sample we consider is made of cetyltrimethylam- 
monium bromide (CTAB) at 0.3 M with sodium ni- 
trate (NaN0 3 ) at 0.405 M. The temperature is fixed at 
T=28°C. This system is Maxwellian in the linear regime, 
with a single relaxation time tr = 0.23 ± 0.02 s and a 
plateau modulus Go = 238 i 5 Pa while, under simple 
shear flow, it is well-known to exhibit shear banding asso- 
ciated with an instability of the interface and Taylor-like 
vortices [1, [13, E3 ■ Our experiments are performed in 
a cylindrical Couette device with inner rotating cylinder 
(height H = 40 mm, inner radius Ri — 13.33 mm, gap 
e = 1.13 mm) Q adapted to a stress-controlled rheome- 
ter (Physica MCR301) used in strain-controlled mode. 
Fig.Q]displays the overall rheological behavior of the sam- 
ple for Wi between and 34 together with simultane- 
ous observations in the velocity gradient-vorticity (r, z) 
plane summarizing the main flow states along the flow 
curve (inset). For Wi < Wii, the gap of the Couette 
cell appears homogeneous (inset-a) and the correspond- 
ing branch of the flow curve is slightly shear-thinning. 
Between Wii — 1.0 ± 0.1 and Wih = 26 ± 1, the shear 
stress presents a plateau and the sample splits into two 
shear bands of differing optical properties separated by 
an interface undulating along the vorticity direction. The 
flow is organized in Taylor-like vortices, mainly localized 
in the high shear rate band (inset-b) (8|- Note that the 
value of Wih is given from the crossover in the flow curve 



between the stress plateau and the high shear branch 
and is in good agreement with the Weissenberg num- 
ber for which the high shear rate band fills the whole 
gap. The standard deviation is computed from statisti- 
cal measurements (see Ref. [HI for more details). Let 
us now focus on the high shear rate branch of the flow 
curve (Wi > Wih), where the induced high shear rate 
phase fills the entire gap. Two different regimes can be 
distinguished : (1) For Wi between Wih and a critical 
value noted Wi c — 30 ± 0.5, the flow curve is strongly 
shear-thinning and vortices seem to disappear (inset-c). 
All the methods used in ref. Q to identify vortices in the 
banding regime have failed. The flow in this regime seems 
purely ID. This behavior is compatible with the observa- 
tions of Salmon and coworkers, who demonstrated using 
velocity measurements on a similar system that the in- 
duced structures were highly non-Newtonian. (2) For 
Wi > Wi c , the flow curve presents a drastic change of 
variation corresponding to an apparent shear-thickening 
behavior. This increase of the flow resistance is due to 
the new structure of the flow [l3|], and differs from the 
nucleation of a new structure in the material itself fl9| . 
Indeed, as seen in Fig. [T] (insct-d) and in the supple- 
mental movie [20| . in this regime, the flow appears very 
disorganized in space and time. This behavior is strongly 
reminiscent of the state of elastic turbulence highlighted 
by Groisman et at in dilute polymer solutions flowing in 
different flow geometries [l]| . According to Ref. [l]| , the 
main features of the elastic turbulence are the following : 
fluid motion excited in a broad range of spatial and tem- 
poral scales, and significant increase of the flow resistance 
and rate of mixing. Furthermore, the transition to the 
turbulent regime is found to be hysteretic and strongly 
subcritical. Using mechanical measurements and direct 
visualizations of the sample in the (r, z) plane, we care- 
fully analyse the irregular flow we observe. 
Fig. [2^ shows the rheological behavior of the high shear 
rate branch of the flow curve during increasing and de- 
creasing shear rate sweeps. When the shear rate is 
changed quickly, the transition exhibits a pronounced 
hysteresis. Slower shear rate sweep leads to a re- 
duced hysteresis loop, as expected for subcritical tran- 
sitions [21| . Another strong indication of the subcritical 
character of the transition is given in Fig. [2Jd where, for 
imposed Weissenberg numbers in the transition range, 
the system exhibits bi-stability. As observed for polymer 
solutions, the transition towards the apparent disordered 
state is characterized by a sudden stress jump. The visual 
impression of spatiotcmporal disorder in Fig. [T]inset-d is 
confirmed by a more careful analysis. Fig. [3]Ja) shows 
average Fourier spectra of the intensity profiles along the 
vorticity direction, which exhibit power law decay over 
two decades in the wavevector domain. The power spec- 
trum in wavevector seems to follow a power law up to 
kd ~ 100 mm -1 . This suggests that for lengths below 10 
/xm, the dissipation process is different, leading to a faster 
relaxation of the fluctuations. In the case of elastic tur- 
bulence, kd is expected to be linked to the stress diffusion 
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FIG. 2: (a) Close-up on the high shear rate branch of the 
flow curve (Wi > Win) obtained for increasing (circles) and 
decreasing (squares) shear rate, for two different sampling 
times : 2 s/point (close symbols) and 1 min/point (open sym- 
bols). a u and ad are the apparent up and down boundaries 
of the stress jump at Wi c . (b) Shear stress evolution on long 
time (between 1 min and 1 hour), for different Wi. From 
top to bottom: Wi = 32.8 (fully turbulent), Wi = 31.2 (bi- 
stable), Wi = 30.1 (fully laminar). 



coefficient by fcj" ~ V S#tr [22(. When constitutive 
models are used to describe the dynamics of surfactant 
solutions, 2> is typically the coefficient of non-local diffu- 
sive terms . Hence from the value of kd we observe here 
in the turbulent regime, we can compute a stress correla- 
tion length of the order of 10 /xm. This estimation agrees 
with the stress correlation length (4 to 8 /xm) obtained 
from completely different experiments in Ref. fl8ll23l. |24| . 
Similarly, the intensity reflected by the sample at one 
point over time gives us some information about tem- 
poral disorder in the system. Fig. E^b) shows average 
Fourier spectra of the intensity at different points in the 
gap over time. The spectra exhibit power law decay over 
two decades in the frequency domain. Furthermore, in 
order to get information on the fluctuations at the global 
scale, we have studied the time series of the macroscopic 
shear rate measured by the rheometer when imposing 
constant stresses. The inset of Fig. [1] gives the evolu- 
tion of the amplitude of the shear rate fluctuations as a 
function of the shear stress. Below the transition to tur- 
bulence, fluctuations are small (AWi=l), less than 3% 
of the average value. In the transition range fluctuations 
are very high (AW?=11), up to 30% of the average value. 
In the fully turbulent regime, i.e. <j/Gq > 1, fluctuations 
are as high as 15% of the average value (AWi=5). Av- 
erage Fourier spectra of the fluctuations in this regime 
exhibit power law decay over almost two decades in the 
frequency domain, with an exponent j3 of —3.5. This 
value of the exponent satisfies the criterion derived by 
Fouxon et al. (/3 < —3) who have rationalized the spec- 
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FIG. 3: Average Fourier spectra of the intensity reflected by 
the sample in the velocity gradient-vorticity plane for two Wi 
larger than Wi c . Circles: Wi — 31.5. Squares: Wi — 33. 

(a) Spectrum obtained on the wavevector domain by com- 
puting Fourier transforms of the intensity reflected along the 
vertical direction at a midle point in the gap and then by 
averaging for different times. The dashed line is a k~ 2 fit. 

(b) Spectrum obtained on the frequency domain by comput- 
ing Fourier transforms of the intensity reflected overtime at 
single points in the gap and then by averaging for different 
points. A fast camera (Phantom V9) has been used to reach 
the frequency range between 10 Hz and 700 Hz. The dashed 
line is a f~ 1A fit. 



tra of turbulence in polymer solutions [22| . 
Hence, all the features highlighted for Wi > Wi c indi- 
cate that the dynamics in this regime are likely due to 
the same mechanisms driving the elastic turbulence in 
polymer solutions. Note however that the transitional 
pathway between the (ID) laminar and the turbulent 
regime in polymer solutions appears different. In par- 
ticular, when decreasing Wi, we do not observe solitary 
vortex pairs [25j , the relaxation of turbulence being ho- 
mogeneous. Differences might arise from the existence, 
for micelles, of two distinct relaxation mechanisms un- 
derneath tr (r r and T),). 

Let us summarize the complete scenario we observe : 
in the shear banding regime, the flow is unstable with 
respect to perturbations along the vorticity direction 
and Taylor-like vortices develop in the high shear rate 
band. The stability of the flow seems to be recovered 
for Wih < Wi < Wi c , namely when the high shear rate 
phase fully fills the gap, up to a critical Weissenberg num- 
ber Wi c above which the system undergoes a transition 
towards a turbulent state. Taking into account the clas- 
sical succession of instabilites expected with increasing 
Wi [HI, EH , from a purely ID flow to a non trivial coher- 
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FIG. 4: Average Fourier spectra of the macroscopic shear rate 
obtained by imposing three constant stresses in the turbulent 
regime. The total duration of each time serie is 5000 s and 
the sampling time is 0.1 s. Triangles : a/Go ~ 1. Circles : 
a/Go — 1.03. Squares : a/Go = 1.05. The dashed line 
is a /~ 3 ' 5 fit. (inset) Average amplitude of the fluctuations 
of shear rate (i.e. Wi) for different imposed stresses before, 
"during", and after the transition to turbulence. Stress is 
measured in unit of the stress jump in the flow curve at Wi c . 
For stresses between aa and a u (i.e. < ^~"„ d < l)j the 
system is bi-stable and exhibits large fluctuations. Spectra 
are computed in the fully turbulent range, for stresses larger 
than a u . 



ent flow and finally to a turbulent flow, the stability of the 
high shear rate phase in the range Win < Wi < Wi c is 
intriguing. One explanation could be linked with the fact 
that instability thresholds strongly depend on boundary 
conditions. From the point of view of the high shear rate 
phase, the boundary conditions change with increasing 
Wi. Indeed, during the shear banding regime, the high 
shear rate phase is confined between a rigid wall (the 
inner rotating cylinder) and the viscous band that acts 



as a "soft boundary". In contrast, for Wi > Wih, the 
high shear rate phase is in direct contact with the two 
rigid walls of the Couette cell. In the same manner as 
many other instability examples, we can expect the in- 
stability threshold to be lower for soft boundaries than 
for rigid boundaries (26j . A more precise picture follow- 
ing this trend could be given by an analysis focused on 
the high shear rate phase and considering the viscous 
phase as an elastic boundary [27} • Strictly, the interme- 
diate stable regime may also suggest that instability in 
the plateau regime is due to the interfacial mechanism, 
and that the bulk instability of the high shear rate branch 
for Wi > Wi c is a classical elastic turbulence. The sta- 
bility analyses that uncovered the interfacial mechanism 
have been performed in planar geometries. Very recently 
S. Fielding extended her original calculation to Taylor- 
Couette flow to test the effect of curvature on the inter- 
facial mechanism (28j . Using a unique framework, she 
suggests that both interfacial and bulk elastic mecha- 
nisms could be observed for different values of stream- 
lines curvature and normal stresses. This first theoretical 
work dealing with shear banding and elastic instabilities 
together predicts that in a curved geometry similar to 
ours, the interfacial mechanism is less favorable. In con- 
clusion, our findings establish clearly that elastic insta- 
bilities and shear banding are a lot more intertwined than 
presupposed. These findings, together with the most re- 
cent theoretical studies, now rationalize 3D disturbances 
during shear banding in flows with curved streamlines as 
being driven by elastic instability. 
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